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Abstract 


In this paper, the variational iteration method which proposed by Ji-Huan 
He is applied to solve both linear and nonlinear Schrodinger equations. The 
main property of the method is in its flexibility and ability to solve linear 
and nonlinear equations accurately and conveniently. In this method, gen- 
eral Lagrange multipliers are introduced to construct correction functionals 
to the problems. The multipliers in the functionals can be identified opti- 
mally via the variational theory. Numerical results show that this method 
can readily be implemented with excellent accuracy to linear and nonlinear 
Schrodinger equations. This technique can be extended to higher dimensions 


linear and nonlinear Schrodinger equations without a serious difficulties. 
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1 Introduction 


The variational iteration method (VIM) was first proposed by Ji-Huan He in 1998 
[6,7] and systematically illustrated in 1999 [11]. Since then, it has been success- 
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fully applied to various engineering problems [15,16]. This method is employed 
in [18] to solve the Klein-Gordon equation which is the relativistic version of the 
Schrodinger equation, which is used to describe spinless particles. Application of 
He’s variational iteration technique to an inverse parabolic problem is described 
in [4]. In [2] the VIM is employed to solve the time dependent reaction-diffusion 
equation which has special importance in engineering and sciences and consti- 
tutes a good model for many systems in various fields. This technique is also 
employed in [5] to solve the Fokker-Planck equation and in [3] to solve a bio- 
logical population model. For more application of the method, the interested 
reader is referred to [1,17,19,21]. The VIM [11,12] is a powerful tool to search 
for approximate solutions of linear and nonlinear equations without requirement 
of linearization or perturbation. Another important advantage is that the VIM 
is capable of greatly reducing the size of calculation while still maintaining high 
accuracy of the numerical solution. Moreover, the power of the method gives a 
wider applicability in handling a huge number of analytical and numerical ap- 
plications. The convergence of He’variational iterative method is investigated 
in [20]. Here, we apply VIM to one and two dimensional linear and nonlinear 
Schrodinger equations. This paper is organized as follows: In Section 2, we in- 
troduce the model of the problems. In Section 3, first we describe VIM method 
and then we apply VIM in a direct manner to establish exact solutions for linear 
and nonlinear Schrodinger equations. In Section 4, we describe the numerical 
solution of linear and nonlinear Schrodinger equations to show the power of the 


method in a unified manner without requiring any additional restriction. 


2 The model of the problem 


In this paper, the linear Schrodinger equation is considered as follows: 


a (x,t) = 90, w(x, 0) = f(z), ceER, t > 0, aed, (1) 
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and we consider the nonlinear Schrodinger equation of the form 


OW 1 Ay Pap 2 2 
ix 3 (age (8) + Gal%b) + Bly, XER', £20, (2) 


where X = (x,y) , |4|? = a, and £ is a real constant. 


(X, t) ee 


3 Basic ideas of He’s variational iteration method 


In this section, the application of the VIM is discussed for linear and nonlinear 


Schrodinger equations. Considering the following general differential equation: 

Lnp(a, t) + Roa, t) + Nep(a, t) = g(2,t), (3) 
where L is a first order partial differential operator, # is a linear operator, N is 
a nonlinear operator and g(z,t) is a known analytical function. According to the 
VIM[8-10], we can construct the following correction functional: 


t ~ ~ 
Yn+i (2, t) = in(e,t)+ f ALin(f) + Rin (§) +NYn(S)—9(8))d& nn 2 0, (A) 


0 
where 4 is a general Lagrange multiplier [14], which should be identified optimally 


via the variational theory [14], the subscript n denotes the nth approximation, 
and wn is considered as a restricted variation [6,7,11] and [13] ic dy, =0. We 
first determine the Lagrange multiplier A that will be identified optimally via 
integration by parts. The successive approximations ~p4+41(x,t) ,n > 0 of the 
solution w(z,t) will be readily obtained using the derived Lagrange multiplier 
and by using any selective function wo. The initial values w~(z,0) and 7;,(z, 0) 
are usually used for selecting the zeroth approximation wo . With 4 determined, 
several approximation ~j;(z,0), j > 0 follow immediately. Consequently, the 


exact solution may be obtained by using (see [20]) 

= lim Yn. (5) 
According to the VIM, we consider linear Schrodinger equation (1) in the follow- 
ing form[8—10]: 


On 


PH (0, 8) +95 (a Ode. (6) 


og 


t 
Png (t,t) = Ynle,t) + [ AE) 
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To find the optimal value of A, we have 


t Qa) 
Fran ast) = durn(ast) +5 f XE FE (e.8) +4 FM. ))ds =0, (7 
or - ; os - 
Firnsi(e.t) = Siba(v.t) +5 f° E(w. €))ag = 0. () 


which follows 


t 
Hnsi(tat) = bbu(e,t)(1 +) 5 f M@vale.g)dé=0, (9) 
The following stationary conditions 


Lt Kay S 0 (10) 


follow immediately. This in turn gives 


A(é) = -1. (12) 


Substituting this value of the Lagrange multiplier \ = —1 into the functional (6) 
gives in the following iteration formula 


On 


“Pn (6) +4 EE (0). (18) 


og 


t 
saat) = wales) — f( 
Similarly, we obtain the correction functional for (2). Hence we have 


Orn 1 Pun un 


inl) = onlXst)+ [xO Beag + (Se w9 + Hx.) 
Bp p)dé. (14) 

The stationary conditions are of the following form 
1+ éX(t) =0, (15) 
N(é) = 0, (16) 


and so we have 
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Substituting this value of the Lagrange multiplier A = 7 into the functional (14) 


gives the following iteration formula 


t 6 Ms 2 e bn 
dalXt) = dalXst) +6 f OPES + (GEO + Fe) 


Bp pdé. (18) 


Here, we will use this method to solve linear and nonlinear Schrodinger equations 


to establish exact solutions for these equations. 


4 Examples 


To illustrate the solution procedure and show the ability of the method, some 
examples are provided. 


Example 4.1 Consider the following linear Schrodinger equation : 


a e 
Sola, t) +o (ut) = 0, (19) 
wo(x) = sinh 2a. (20) 


Using (13), we obtain the following successive approximations: 


Yi(z,t) = (1— 4:t) sinh2z, 
—4it)? 
ities (=. tray - JF Atos 
—4it)?  (—4it)? —4it)” 
ei = @aaees x been i JE he hl OE ieee 
Consequently, the exact solution is 
Wenact(£, t) = e * sinh2z. (21) 
Example 4.2 Consider the nonlinear Schrodinger equation 
OW Pap 2 
iP (0,1) + 3 (w, 1) +2|bPy = 0, (22) 


po(x) = pt (23) 
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Using (18), we obtain the following successive approximations: 


di(@,t) = (1+it)e™, 
+ 2 : 
volet) = +) + Dyer, 
a ae Ci ee 
vs(et) = + G+ Gy Bye, 
4) 2 3 4) 4 
Yn(z,t) = (1+ (it) 4 “0 “ “ t ... 4 “ je™ 
Consequently, the exact solution is 
enact (yt) eh), (24) 
Example 4.3 Consider the nonlinear Schrodinger equation 
OW 1 A*y O2yp 
Pag OO Gh aan AP a) + 2|p|?y =0, (25) 
Po(x,y) = bY, (26) 
Using (9), we obtain the following successive approximations: 
Wi(X,t) = (1+ itye(@*), 
. 2 . 
qo(X,t) = (14 (it) CEP eile, 
. 2 . 3 . 
wiX) = C+ n+ Se 4 Seer, 
ab) GES». Gaby? ee 
un(X,t) = (1+ (it)4 “ “ “ fea ced “ jelly) 


Consequently, the exact solution is 


Weaact(X, t) = ellttety) 


Conclusions 


In this paper, He’s variational iteration method has been successfully applied to 


find the solution of the linear and nonlinear Schrodinger equations. The main 
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advantage of the method is the fact that it provides an analytical approximation, 
in many cases an exact solution, in a rapidly convergent sequence with elegantly 
computed term. Analytical solutions enable researchers to study the effect of 
different variables or parameters on the function under study easily. A clear 
conclusion can be drawn from the numerical results which VIM provides with 
highly accurate numerical solution without spatial discretizations for linear and 


nonlinear Schrodinger equations. 
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